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In 1955 Linde Air Products introduced a plasma arc torch for 
metal cutting applications, and by 1965 Linde had developed an 
automatic plasma arc welding facility for Westinghouse Electric 
Corporation to fabricate 120 in. diameter, 3/8 in. thick D6AC steel 
rocket motor case for the Titan III-C booster assembly (Miller and 
Filipski, 1966; Privoznik and Miller, 1966). The plasma arc welding 
process was reported to have halved the welding time required. Then 
in 1965 Thermal Dynamics Corporation reported the use of direct 
current reversed polarity (workpiece negative) plasma arc to join 1/4 
in. thick aluminum plate (Cooper el al., 1965). At the end of the 
decade of the 1960's, Van Cleave at the Boeing Company began his 
efforts to combine the plasma arc process with a variable polarity 
feature wherein the electrode polarity was periodically reversed. 
Alternating electrode potential for aluminum welding had been 
investigated as early as 1947 (Herbst, 1947) . Difficulties with 
welding power supplies in this application were evident early-on and 
when variable polarity plasma arc (VPPA) welding was used in the U.S. 
Army Roland Missile Production Program, development problems such as 
arc pressure pulsation were noted (Nunes et al., 1981; 1983) . As a 
result of Van Cleave's promising work at Boeing, a VPPA research and 
development project was initiated in 1978 at the NASA Marshall Space 
Flight Center, intended to determine the potential for replacing the 
GTAW system used in the fabrication of the Space Shuttle External Tank 
(ET) . Certain improvements were made to the original equipment and 
process, VPPA welding finally exceeded the expectations (Nunes et 
al., 1984a; 1984b). 
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Ths Space Shuttle ET is the largest known "drop tank" carrying 
140,000 gallons of liquid oxygen and 380,000 gallons of liquid 
hydrogen, with the dimensions 28.6 feet in diameter and 154 feet in 
length. From the initiation of ET production, the conventional GTAW 
system operating in the direct current electrode negative mode was 
used. In fabricating the large Saturn lunar rocket first stage with 
the GTAW system, weld porosity and inclusions had been ever-present 
P^°klsms throughout the entire program, and the decision to use 100% 
radiographic inspection was made at the outset. Weld porosity and 
inclusions had to be systematically ground out and weld repaired which 
made the replacement of GTAW system by VPPA welding process 
advantageous . 

Significant advantages of the VPPA process include faster 
welding, fewer repairs, less joint preparation, reduced weldment 
distortion, and absence of porosity. 

In this study, flow profile and power distribution of argon arc 
plasma through the entire VPPA welding process has been analyzed. 
This study includes the system analysis of following sections: (l) 

and main arc power distribution and the total power input of 
the VPPA welding system; (2) Plasma arc power loss through nozzle 
which accelerate arc plasma gas and transform potential energy to 
kinetic energy which form the plasma jet at the outlet of nozzle; (3) 
Radiative and convective heat loss of plasma jet column from the 
outlet of nozzle to the plasma jet approaching workpiece; (4) Flow 
criteria of partial penetration and full penetration plasma arc jet 
through workpiece; and (5) Determination of crown width, crown 
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height, root width and root height for horizontal and vertical 
welding. 

II. Variable Polarity Electric Power System 
The significant current carriers in a welding arc are electrons 
and positive ions. Electrons carry the bulk of the current, moving 
rapidly from negative cathode to positive anode. The positive ions 
more slowly through the interelectrode space. The 
drift rate results in an unsymmetrical heating at the 
ends of the welding arc in a fixed polarity arrangement. The cathode 
receives less heat and the anode more heat. The "straight polarity" 
mode of operation entails a negative electrode (cathode) and a 
positive workpiece (anode) . Where the primary object of the weld 
process is to deliver the maximum heat to the workpiece with minimum 
deterioration of the electrode, straight polarity is used. 

Reverse polarity", which is positive electrode and negative 
workpiece, has the advantage that the workpiece is subjucted to a 
cleaning process (cathode cleaning) by the impingement of heavy 
(compared to electron) positive ions on the workpiece surface. In 
the case of plasma arc welding reverse polarity action appears to 
condition the surface of the aluminum alloy so that the molten metal 
flows easily and controllably under the arc. It is conjectured that 
this fluid control is accomplished through breaking up of surface 
oxide films. Cutoff of reverse polarity during VPPA welding 
transforms a weld metal flow which closes smoothly and soundly behind 
the keyhold to an irregular, intermittent, globular flow leaving a 
rough, lumpy bead sunken below the parent metal surface and protruding 
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jaggedly from the root of the weld. But continuous reverse polarity 
is not necessary to provide adequate cathodic cleaning action. 

The variable polarity square-wave with unequal straight and 
reverse polarity time intervals offered a combination of the high 
heating capability of straight polarity with the cleaning feature of 
reverse polarity. Adequate cleaning is obtained by incorporating a 
relatively short (one tenth to one fifth the duration of the straight 

polarity current) pulse of reverse polarity current into welding 
current waveform. 

Figure 1 show the VPPA electric power system. During the 
straight polarity mode of operation, it employs a negative electrode 
(cathode) and a positive workpiece (anode) , while it employs a 
positive electrode and a negative workpiece during the reverse 
polarity mode of operation. At left of Figure 1 illsustrates the 
straight polarity mode, and the reverse polarity mode is shown on the 
right-hand- side of the figure* 

Let i> x , <p Zl and <t> 2 denote the workfunctions for tungsten 
electrode, nozzle and workpiece, respectively. v a and V c are anode 
and cathode drops, respectively; while ( f + )i and ( 0l are the pilot 
arc voltages in straight polarity and reverse polarity, respectively, 
u shows the voltage drop per unit length in the arc column. One can 

derive the pilot arc voltages in the straight polarity and reverse 
polarity as follows: 

(A) Straight Polarity 

U E = U F - V a - = -V a - <t> 3 

U c = U E - u(L x + d 3 ) = -v a - , i 2 - u ( L x + d 3 ) 




U B = U E -U(d 3 + d 2 ) = -V a - * 3 - u(d 3 + dj) 

U D * U c + v a + 4> z = -v a - i> 2 - u(L t + d 3 ) + v a + <p 2 

- *2 - *3 - U ( L 1 + d 3 ) 

U A = U B - V C + 5*1 = “ V a “ *3 “ u ( d 3 + d 2) ~ V C + * i 

(f+)l = U D ~ U a = 0 2 “ 0 3 ~ u (Lj + d 3 ) + V a + ^ 3 

+ u(d 3 + d 2 ) + v c - <t> 1 

= 5* 2 “ *1 + V a + V C + u ( d 2 - L : ) 

(B) Reverse Polarity 

Voltage difference between point E and E', shown in Figure 1 is 
affected by the cathode section of workpiece. As the aluminum is 
not an ideal thermionic material, the voltage drop of workpiece could 
be greater than that of the ideal thermionic metal, such as tungsten. 
Assuming U EE ' = a, pilot arc power supply in reverse polarity can be 
evaluated as follows: 


U E • - U f • + 

V C - 03 - V c - 

^ 3 



U c ’ = U E • + 

a + u (d 3 + Lj ) 

= V c - i> 2 

+ 

U (d 3 + Lj ) + a 

U B - = U E • + 

a + u(d 3 + d 2 ) 

- V C - *3 

+ 

u(d 3 + d 2 ) + a 

u 0 * = U c • + 

V a + *2 - V c - 

s6 3 + u(d 3 

+ 

L| ) + a + V a + 

U»* = U B * + 

V a + *i « V c - 

* 3 + u ( d 3 

+ 

d 2 ) + a + V a 4* 

(€.) 1 = U D • 

- Ua‘ 




= u c - 

0 3 U ( d 3 Lj 

) + a 4- V a 

+ 

0 2 - V c + <t>2 

~ u(d 3 + d 2 ) -a- 

V a - *1 



I 

(V 

II 

*i - u(d 2 -Li) 





Reasonable values obtained from the literature and from limited 
experimentation at NASA Marshall Space Flight Center yield the 
following values specified: <j> x = 2.6V, </> 2 = 4.5V, <f> 2 = 4.2V, V c = 
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6.5V, V = 7.0V, u= 35.5V/in. Typical geometric parameters are : L : 

= 0.133 in, L 2 = 0.375 in, d 2 = 0.178 in, and d 3 = 0.22 in. We can 
calculate the pilot arc voltages in the straight polarity and the 
reverse polarity as follows: 

( € + ) 1 = i> 2 ~ 01 + V a + V C + u ( d 2 " L l) 

= 4.5 - 2.6 + 7.0 + 6.5 + 35.5(0.178 - 0.133) 

= 17 V 

( O 1 = *2 " ~ u ^ d 2 “ L l) 

= 4.5 - 2.6 - 35.5(0.178 - 0.133) 

= 0.3 V 

III. Total Power Input of the VPPA Welding System 

The pilot arc current, superimposed upon the main arc current 
between the electrode and the plasma nozzle, is supplied by a seperate 
power supply located in the plasma control console. The pilot arc 
polarity does not alternate. The nozzle is held positive with 
respect to the electrode so that the electron flow is always away from 
the electrode toward the nozzle. The main arc is started by 
application of high frequency alternating voltage between the 
electrode and the plasma nozzle so as to first establish the pilot arc 
from which the main arc proceeds. The high frequency current is 
controlled by the plasma control console, but is generated by the 
spark gap oscillator located in the hot block. 

The total electric power input of VPPA welding system includes 
the pilot arc between electrode and nozzle, and also the main power arc 
between electrode and workpiece during the straight polarity and the 
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reverse polarity modes. Following nomenclatures are given: 

I x = pilot arc current 
I 3 = main arc current 
= pilot arc voltage 
e 3 = main arc voltage 

(I + ) 3 = main arc current in straight polarity 
(I _) 3 = main arc current in reverse polarity 
t + = time period of straight polarity 
t_ = time period of reverse polarity 

The total electric power input of the VPPA welding system is 
given in the expression 

P = I x e : + I 3 € 3 

( f +)l^ + + (€_)jt_ (i+) 3 ( f +) 3 i-+ “*■ (^— ) 3 ( ^_) 3 *t-_ 

= Ij + 

t + + t_ t + + t_ 

Figure 2 shows the power input system and the dimension of the 
geometry. In the present VPPA welding system operated in the 
NASA/MS FC , time periods for the straight polarity and the reverse 
polarity are 19 sec and 4 sec, respcectively . Pilot arc currents for 
both straight and reverse polarities are 25 amp; while main arc 
currents for straight and reverse polarities are 195 and 200 amp, 
respectively. The pilot arc voltages for the straight and the 
reverse polarities are 17 and 0.3 V, respectively; while the main arc 
voltages for the straight and the reverse polarities are 32.5 and 40 V, 
respectively. Thus, we can calculate the total power input of the 
VPPA welding system operated in the NASA/MSFC as follows: 
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(17) (19) + (0.3) (4) 


P = (25) 

(19) + (4) 

(195) (32.5) (19) + (200) (40) (4) 

+ 

(19) + (4) 

= 352.5 + 6626.63 = 6979.13 W 

The volumetric flow rate of argon gas to be used for plasma arc 
jet stream is regulated. In the present VPPA welding system operated 
in the NASA/MSFC , volumetric argon gas flow rate is 7 . 5 ft 3 /hr at 25°C. 
It can be converted to SI unit as follows: 

Q = (7.5 ft 3 /hr) (1/3600 hr/s) (1/35.3147 m 3 /ft 3 ) • 

(10 6 cm 3 /m 3 ) 

= 58.99 cm 3 /s 

Argon gas density at 1 atmospheric pressure and 25 °C can be computed 
from the equation of state 

R 

P = p — T 
M 

where p denotes pressure; p, the density; R (=82.05 cm 3 * atm/mole • °K) , 
the universal gas constant; M, the molecular weight (39.948 g/mo le for 
argon gas) ; and T, the temperature. The argon gas density at 1 
atmospheric pressure and the temperature of 25 °C can be computed as 

(82.05 cm 3 * atm/mole* °K) 

(1 atm) = p — ( 298 °K) 

(39.948 g/mole) 

which leads 

p = 1.6249 x 10 -3 g/cm 
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Mass flow rate of argon gas, thus, can be obtained as 

• • 

m = pQ = (1.6249 x 10‘ 3 g/cm 3 ) (58.99 cm 3 /s) 

= 9.585 x 1CT 2 g/s 

Shielding gas is needed for the isolation of welding system from 
the surrounding environment, and for the cooling of workpiece from the 
high temperature injection of plasma arc. Shielding gas is an 
essential factor in achieving quality welds in any gas shielded arc 
welding process. The problem of impurities, when manifest, is 
normally resolved by eliminating the source of the contamination. 
However, if the gas supply itself is contaminated serious schedule 
interruptions can result, particularly where large manifolded gas 
distribution systems, are employed. 

Gas contamination affects VPPA welding in several ways. Inside 
the torch, oxygen from water vapor decomposed in the arc can produce a 
visible oxide film stain on the tungsten electrode. The plasma jet 
may sputter and exhibit a green coloration from copper eroded from the 
nozzle under erratic flow conditions. The weld itself may exhibit a 
pocked- looking surface, presumed to be fine porosity due to hydrogen 
gas from decomposed water vapor. Under more severe conditions, the 
torch may leave behind it a train of irregularities and holes in the 
weld bead. 

IV. Plasma Arc Power Loss Through Nozzle 

Plasma arc power distribution in the whole variable polarity 
plasma arc welding system can be divided into the following four parts 
for the convenience of discussion: (1) plasma arc power before 

entering the nozzle (P a ) ; (2) plasma arc power after extracting from 
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the nozzle (P b ) ; (3) plasma arc power before entering the keyhole of 
workpiece (P c ) ; and (4) plasma arc power at the outlet end of the 
keyhole of workpiece (P d ). Figure 3 shows the plasma arc power 
distribution in the whole VPPA welding system. 

In this section, plasma arc power loss through the nozzle will be 
discussed. Figure 4 shows the electric current and electric voltage 
drop distribution of the VPPA welding system. In addition to the 
nomenclatures given in Section III, there are some more terms, shown 
in Figure 4, are given as follows: 

Ae + = voltage drop of electrode in straight polarity 

Ae_ = voltage drop of electrode in reverse polarity 

€ + ' = voltage drop of nozzle in straight polarity 

e -' = voltage drop of nozzle in reverse polarity 

Estimation of plasma arc power distribution is based on the 
linear power distribution assumption of argon plasma gases between 
the point of entering the nozzle and the point of entering the keyhole 
of workpiece. Plasma arc power before entering the nozzle can be 
computed from the following equation: 


P a = I. 


(A<f + ) : t + + (AC_)it_ (I + ) 3 U* + ) 3 t + + (I.) 3 (AO 3 t_ 


+ 


t + + t_ 

( I + ) 3 (f + , ) 3 t + + d-) 3 («-')t. 


t + + t _ 
( d 2 ~ Ll) 


t + + t. 


ld 2 + d 3 / 


= (25) 


(3.9) (19) +(9.6) (4) (195) (3.9) (19) +(200) (9.6) (4) 


( 19 ) + ( 4 ) ( 19 ) + ( 4 ) 

(195) (14. 1) (19) +(200) (27.4) (4) (0.178-0.133 


( 19 ) + ( 4 ) 


(0.178+0.22 
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= 122.28 + 962.15 + ( 3 224 . 4 ) ( 0 . 113 ) = 1448.79 W (4-1) 

o 

Plasma arc power loss through the nozzle (q R ) can be 
computed from the summation of heat convection, h # A(T plasma - T w ) ; 
free electron loss, (I 1 /e) [ (3/2) kT*] ; and power loss due to 
difference in anode potential and plasma column drop of argon gas at 
nozzle, X i ^ i , namely 


<3a 


h R A(Tpi asma T h ) + 


(3 

-kT* 

2 


+ I 


l^i 


(4-2) 


where h fl stands the film coefficient; A, the circumferential area of 
the nozzle; T plasma , the plasma gas temperature; T w , the wall 
temperature of the nozzle; e, electric charge of electron; k (=1.38 x 
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erg/mole • °K) , the Boltzmann constant; T* , the average 


temperature between plasma gas and the surface of nozzle; and $ % , the 
average difference in anode potential and plasma column drop of argon 
gas . 


(A) Calculation of Film Coefficient h fl 

Heat transfer equation for plasma flow in a tube from Hsu and 
Rubinsky (1987) gives 

N u = 0.2233 (X + )-°- 7455 + 3.66 (4-3) 

where 

x 1 

X + (4-4) 

R ( R ePr)m 

hD 

N u = — = Nussel Number (4-5) 

k 

uD 

R e = — = Reynolds Number (4-6) 

u 
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c p u 

P r = = Prandtl Number 


(4-7) 


Here, x denotes the distance in axial coordinate; R, the radius of 
tube; h, the convective heat transfer coefficient (or film 
coefficient); D (=2R) , the diameter of the tube; k, the heat 
conduction coefficient; u, the flow velocity of plasma gas; u, the 
viscous coefficient; and C p , the constant pressure specific heat. 
Substituting Equations (4-4) and (4-5) in Equation (4-3), we have 


h = — (0.1116) 

R l 


r R R e P r 


0 . 7 4 5 S 


+ 1.83 


(4-8) 


Define h as the average heat transfer coefficient along the axial 
distance of nozzle. After substituting Equations (4-6) and (4-7) in 
Equation (4-8) , it leads 


(0.4385) 


0 . 7 4 5 S 


+ 1.83 


(4-9) 


(B) Determination of Plasma Temperature at the Nozzle 

Entrance of Argon Plasma Gas 

a 

Argon Plasma gas enthalpy (h) p at the nozzle entrance can be 
obtained as 


(h) P - — - 


1448.79 J/s 
9.585 X 10' 2 g/s 


(0.239 cal/J) 


= 3612.5 cal/g 


(4-10) 


As the enthalpy is a function of temperature, the argon plasma gas 


temperature T, at the nozzle entrance can be determined from the Table 
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shown in Drellishak et al. (1963), namely 
T a = 12600 °K 

(C) Plasma Arc Power Loss Due to Convective Heat Transfer 
In order to determine the convective heat transfer loss, one has 
to understand the temperature profile of flow field. In other words, 
the iteration process has to carry out by assuming the temperature 
profile of the flow field for the calculation of convective heat 
trasfer loss. Temperature profile of the flow field can be checked 
from the calculated heat transfer loss, and vice versa. In this case, 
argon plasma jet temperature at the outlet end of nozzle is 
T b = 8760 °K 

^v®r - age temperature of argon plasma flow through the nozzle is 
T = 1/2(12600 + 8760) = 10680 °K 

Coefficients of argon plasma gas properties (at the temperature 

of 10680 °K) and nozzle (surface temperature is assumed to be T u = 1000 

w 

K) can be determined from Drellishak et al (1963) and Lancaster 
(1986), as follows: 

k = 0.85348 J /m • s ■ °K 
C p = 0.5418 cal/g-°K 
u = 2.8592 x 10 -4 kg/m*s 

Thus, we can calculate Reynolds number, Prandtl number, average film 

coefficient, and convective heat transfer loss as follows: 

* • • 
uD puD mD mD 4m 

U- U uA a ( V 4 ) D 2 irD u 
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(4) (9.585 x 10' 2 g/s)(10 -3 kg/g) 

(it) (0.1562in) (2.54cm/in) (2.8592 x 10' 4 kg/m • s) ( 10" 2 m/cm) 
= 107.6 


k 


( 0 . 54 18cal/ g • 0 K) (4 . 1868J/cal) (10 3 g/kg) ( 2 . 8592x10' 4 kg/m • s) 

(0.85348 J/m-s • °K) 


= 0.76 

k f 

h fl = — (0.4385) 

R l 


( 2mC p ' 0 • 7 4 5 5 


irLk 


+ (1.83) 


0.85348 J/in • s • °K 


(0.1562 


( 2 . 54 ) cm 


( 10 ' 2 m/cm) 


(0.4385) 


(2) (9.585 x 10' 2 g/s) (0.5418 cal/g-°K) 

(TT) (0.133 x 2.54 cm) (0.85348 J/m-s«°K) 

-|0 . 7 4 5 5 > 

(4.1868 J/cal) 

• + (1-83) ■ 

(10" 2 m/cm) 

/ 

= 0.41676 J/cm 2 *s- °K 
Convective Heat Loss = h A A(T-T H ) 

= (0.41676) [7T(0. 1562) (0.133) (2.54) 2 ] (10680 - 1000) 
= 1761.75 W 


(D) Plasma Arc Power Loss Due to Free Electrons 
The vapor temperature on the surface of nozzle can be considered 
to be 3000°K (Metcalfe and Quigley, 1975) . Thus, the average 
temperature between argon plasma gas temperature and vapor 
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temperature on the surface of nozzle is 


T* = 1/2(3000 + 10680) = 6840°K 


Free Electron Loss 



f 3 

- kT* 


\ 2 


(25 coulomb/s) 


(1.602 x 10 


- 1 9 


coulomb/molecule) 


' 3 ' 

.2 


(1.38 x 10 


-16 


erg/mol • # K) (6840°K) 


= 2.21 x 10 8 erg/s = 22.1 W 
(E) Plasma Arc Power Loss Due to Difference in Anode 
Potential and Argon Plasma Gas Column Drop at Nozzle 
Average plasma arc power loss due to the difference in anode 
potential and argon plasma gas column drop at nozzle can be shown as 
Ue + ')t + + 


Ii^i = I x 

t + + t. 

where 


(At/) 

= (e 

+ a ) n 

- 

(6 

♦ > 

n 


(*€_') 

- (e 

- a ) n 

- 

(« 

-) 

n 


(<+)„ 

= (€ 

+ ) c , g 

+ 

( 

e d 

) (d 2 - 

- Li) 

( f + a ) n 

= (€ 

+ ) c , e 

+ 

( 

e + 

)i 


<O n 

= (6 

- ) c , g 

- 

( 

f d 

) (d 2 - 

- Li) 

(«-•>„ 

= (6 

- ) c , e 

+ 

( 

6 _ 

>i 



Here, (f+) C/ g denotes argon plasma gas electric potential in the 
neighborhood of electrode (cathode) during the straight polarity; e d , 
column drop of argon plasma gas; (e + ) C/e , electric potential of 
electrode (cathode) during the straight polarity; ( c -)c,g' ar g° n 
plasma gas electric potential in the neighborhood of electrode 


(anode) during the reverse polarity; and ( e_) c e , electric potential 
of electrode (anode) during the reverse polarity. For the VPPA 
welding system operated in the NASA/MSFC, following values are given: 


( e + ) c , g 

= -25.3 V, 

( d = 35.5 V/in 

d 2 - L : 

= 0.045 in, 

( O c , e = -29.2 V 

(Ol = 

17 V 

(O c ,g = 29.7 V 

(O c ,e 

= 39.3 V, 

( f_) l = 0.3 V 

Thus , 



(O n = 

( O c , g + «d ( d 2 

" Li) 


= -25.3 + (35.5) (0.045) = -23.7 V 

( e + a ) n — ( f + ) c , e + ( f + ) i 

= -29.2 + 17 = -12.2 V 

(On = (Oc /g - O d 2 - Li) 

= 29.7 - (35.5) (0.045) = 28.1 V 

( f - a ) n — ( e - ) c , e + ( e - ) i 

= 39.3 + 0.3 = 39.6 V 

(AO) “ (« +a )n ~ («-) n 

= -12.2 + 23.7 = 11.5 V 
(AO) = (e_ a ) n - (e_)„ 

= 39.6 - 28.1 = 1.5V 
Therefore , 

UO)t + + UO)t_ 

Oi - Ii 

t + + t _ 

(11.5) (19) + (1.5) (4) 

= (25)— = 287.5 W 

(19) + (4) 

(F) Plasma Arc Power Loss Through Nozzle 
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q fl = h A A(T-T w ) + — -kT* + 

s 2 i 

= 1761.75 + 22.1 + 287.5 = 2071.35 W 
(G) Percentage of Plasma Arc Power Loss Through Nozzle 
Among the plasma arc power loss in the whole VPPA welding system, 
the percentage of power loss through the nozzle is: 

q, 2071.35 

— = = 29.6% 

P 6979.13 


for the VPPA welding system operated at the NASA/MSFC. 

V. Plasma Arc Power Loss Between the Outlet of Nozzle to 
the Inlet of Keyhole in Workpiece 
Flow pattern of argon plasma gas between the outlet of nozzle to 
the inlet of keyhole in workpiece is a plasma gas column jet stream. 
(A) Plasma Arc Power at the Outlet End of Nozzle 
Estimation of plasma arc power distribution is again based on the 
linear power distribution assumption. This can be calculated as 
follows : 

/ L 

Pb = + 3 ( f , ) 3 - *3 A 

k d 2 d 3 j 

( 0.133 

= 1448.79 + (3224.4) - 2071.35 

10.178 + 0.22 , 

= 454.94 W 

b 

Argon plasma gas enthalpy (h) p at the outlet of nozzle can be 

obtained as 

b P b (454.94 J/s) 

(h) = = (0.239 cal/J) = 1134.4 cal/g 

m (9.585 x 10 ' 2 g/s) 

19 


can be 


Argon plasma gas temperature at the outlet of nozzle, T b , 
determined from table shown in Drellishak et al. (1963), i.e., 

T b = 8750 0 K 

(B) Plasma Arc Power Loss Between the Outlet End of Nozzle 
to The Inlet End of Workpiece 

Following assumptions are made for proceeding the mathematical 
formulation : 

(1) Major power loss in this region consists of radiative heat 
transfer loss and convective heat transfer loss. 

(2) Equally divided plasma flow path into 5 segments between 
the outlet end of nozzle to the inlet end of workpiece 

Calculations procedures for Segment 1 is illustrated. Results of the 
whole flow fields for the five segments are shown in Table 1. 
Computation methods for the rest of four segments are similar to that 
of Segment 1 shown in this report. 

(B-l) Plasma Arc Column Radiative Heat Transfer Loss for 

Segment 1 

Radiative Heat Flux for Segment 1 can be shown in the following 
equation 

Q R1 = fffAj (T plasnU) j - 300 4 ) 

where <r = 5 . 67 x 10" 8 W/m 2 K 4 , is Stef an-Boltzmann constant. Note that 
argon plasma arc temperature is greater than 10 4 K which is several 
orders of magnitude greater than the surrounding medium temperature. 
Therefore, the emmissivity of radiation from the surrounding is 
negligible in the present case. The value of emmissivity can be 
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chosen from Metcalfe and Quigley (1975) which is f = 0.02247. In the 
meanwhile, the surrounding area of the argon plasma arc column is 


Ai = ( TTd j ) 


= [71(0.1562) (2.54) ] 

' 1 

-(0.22) (2.54) 


5) 


5 


= 0.1393 cm 2 

The inlet temperature of the plasma arc column is equal to the outlet 
end of nozzle, 

Tj # e = T b = 87 50 0 K 

The outlet end temperature of Segment 1 plasma arc column T t ^ 0 , has to 
assume. Let us assume T 1/0 = 10350°K. The average temperature of 
Segment 1 is, thus 
_ 1 

T : = — (8750 + 10350) = 9550°K 
2 

The value of C p can be obtained from this temperature, 

Cp (9550° K) = 0.278 cal/g* # K 

Plasma arc column radiative heat loss for Segment 1 is, thus, 

Q R1 = creA 1 (T t 4 - 300 4 ) 

= (5.67 x 10" 8 W/m 2 -K 4 ) (10~ 4 m 2 /cm 2 ) ( 0 . 02 247 ) 

• (0.1393 cm 2 ) ( 9 5 5 0 4 - 300 4 ) 

= 147 . 61 W 

(B-2) Comparison of Plasma Jet Velocity and Shielding 

Gas Velocity 
(a) Plasma Jet Velocity 

At the outlet end of nozzle, argon plasma temperature is T = 
8750 0 K, and density of plasma is p = 5.5355 x 10" 5 g/cm 3 . We can 
determine the cross sectional area of plasma arc column, 
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A = - D 2 = - [ (0.1562) (2.54) ] 2 = 0.12363 cm 2 
4 4 

Argon plasma jet velocity can be calculated from the 
conservation of mass flow rate, m = pu p A. Thus, we have the plasma jet 
velocity 

m (9.585 x 10" 2 g/s)(10“ 2 m/cm) 

u _ - — = 14 0.1 m/s 

P P A (5.5355 x 10" 5 g/ cm 3 ) (o . 12363 cm 2 ) 

(b) Shielding Gas Velocity 

Volumetric flow rate for helium shielding gas of VPPA 

# 

welding system operated by the NASA/MSFC is Q sg = 3 5 ft 3 /hr at 2 5 °C. 
We can convert the present units to cm 3 /s from following relation 


Q 


s g 


(35 ft 3 /h) 


1 

hr/s 

13600 


(35.3147 m 3 /ft 3 ) (10 6 


cm 3 /m 3 ) 


= 275.3 cm 3 /s 

Based on the equation of state, we have 


R 

P = p T 

M 

where M (=4.003 g/mol for helium) denotes molecular weight. One can 

calculate the density of helium from the following relation 

(82.05 cm 3 • atm/mole • °K) 

(1 atm) = p ( 298 0 K) 

(4.003 g/mol) 

Therefore, we have helium density 

p He = 1.637 x 10 -4 g/cm 3 

Mass flow rate of helium is 

m sg = PHeQsg = (275.3 cm 3 /s) (1.637 X 10’ 4 g/cm 3 ) 
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= 4.5071 x 10' 2 g/s 

Cross-sectional flow area of shielding gas is 

7T r 

A sg* = -(°1 2 - D 2 ) = -[ { (0.709) (2. 54) } 2 - { ( 0 . 1562 ) ( 2 . 54 ) } 2 ] 
4 4 

= 2.4235 cm 2 

Again, shielding gas velocity u s can be computed from the 

conservation of mass flow rate, m sg = p He u s A sg *. It leads 

m sg (4.5071 x 10~ 2 g/s)(10~ 2 m/cm) 

u s = = = 1.136 m/s 

pA sg * (1.637 x 10' 4 g/cm 3 ) (2.4235 cm 2 ) 

The ratio of shielding gas velocity to plasma jet velocity is 

u g 1.136 

— = = 8.1 x 10' 3 

u p 140.1 

Therefore, the effect of shielding gas velocity in the flow zone 
calculation of plasma arc jet can be ignored. 

( B-3 ) Plasma Arc Jet Convective Heat Transfer Loss for Segment 1 
Surrounding the high temperature argon plasma jet, there is low 
temperature helium shielding gas. We are not going to encounter with 
the dissociation problem for the shielding gas since helium is 
monoatomic . Based on Ibele (1963) the neutral helium gas begins to 
break up toward ionization at 1 atm in the neighborhood of 13000 to 
15000 °K, and produces helium ions and electrons. In this particular 
case, the consideration of the ionization of helium gas is not 
necessary because the plasma jet temperature is below the temperature 
required for the ionization of helium gas. 

Based on Abramovich (1963) for the arc column convective heat 
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transfer loss for turbulent jet, it shows that the convective heat 
flux is 


Q c 1 = 48A i f>u p C p ^T 

Numerical constant in the equation is resulted from the length scale 
factor of the edge of turbulent jet column. A x is the surrounding 
area of plasma column; u p , the plasma jet velocity; and ^T, the 
temperature difference between plasma jet and surrounding gas. 
Thus, we can calculate plasma arc column convective heat transfer loss 
for segment 1 as follows; 

48 A 1 . 

Qci = 48A lP uC AT = mC 

3600 A c 

48 (0.1393 cm 2 ) 

= (9.585 x 10" 2 g/s) (0.278 cal/g-°K) 

3600 (0.1236 cm 2 ) 

• (4.1868 J/cal) (9550 - 300) = 15.51 Watts 

Here, A : stands cross sectional area of the plasma arc column. 
(B-4) Total Plasma Arc Column Power Loss of Segment 1 
Total plasma arc column power loss of segment 1 is a summation of 
radiative heat transfer and convective heat transfer losses, namely 


• • • 

Ql = Qri + Qci = 147.61 + 15.51 = 163.12 W 

(B-5) Plasma Arc Power at Outlet End of Segment 1 
(d 2 /5) 

Pi — P[j + l3 e ~ Qi 

d 2 + d 3 


= 454.94 + 


(3224 


4) 


0.044 


10.178 + 0.22/ 


163 . 12 
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648.28 W 


Enthalpy of argon plasma jet at the outlet end of Segment 

1 is 

p i (648.28 J/s) 

h i = T~ = (0.239 cal/J) 

m (9.585 x 10 " 2 g/s) 

= 1616.5 cal/g 

Thus, argon plasma jet temperature at the outlet end of Segment 1 can 
be obtained from the table shown in Drellishak et al. (1963), i.e., 
Tj / 0 = 103 70 0 K 

Similarly, the computational methods for physical parameters 
obtained for Segment 1 can be carried out for Segments 2 to 5 . Table 1 

shows the listings of computed physical parameters for Segments 1 to 
5. 

(B- 6 ) Percentage of Plasma Arc Power Loss Between Outlet 
End of Nozzle to Inlet End of Keyhole 

L 

Plasma arc power losses, P i from Segments 1 to 5 are illustrated 

in Table 1 . Percentage of plasma arc power loss between outlet end of 

nozzle to inlet end of keyhole can be determined from the summation of 

power losses from Segments 1 to 5 divided by the total power of the VPPA 

system which is 

L L L L L 
3b _ p i + P 2 + P 3 + P 4 + p 5 1371 

~ — = = 20 % 

P p 6979.13 

VI. Plasma Arc Power Loss Through Keyhole 

Plasma arc jet power at the entrance point of keyhole is 
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P c = P 5 + 


(I + ) 3 (Ae + ') 3 t + + (I_) 3 Ue_') 3 t_ 
t + + t. 


( 6 - 1 ) 


866.5 W + 


(19 + 4 s) 


(195 coulomb/s) (11.2 V)(19 s) + 


+ (200 coulomb/s) (2 . 3 V) (4 s) 


= 8.66.5 W + 1884.17 W = 2750.67 W 
Here, (ac + ') 3 denotes the workfunction and anode drop of workpiece, 
shown in Figure 4 ; and ( a f _ ' ) 3 shows the workfunction and cathode drop 
of workpiece, also shown in Figure 4. For the system operated at 
NASA/MS FC , (if/) 3 is 11.2V and (if.') 3 is 2.3V. 

Enthalpy of argon plasma jet at the entrance point of keyhole can 
be determined by 

P c (2750.67 J/s) 

h c = - — = (0.239 cal/J) = 6858.74 cal/g 

m (9.585 x 10" 2 g/s) 

Thus, argon plasma jet temperature and density at the entrance point 
of keyhole can be obtained from the tables shown in Drellishak et al. 
(1963) , i.e. , 

T c = 14 2 8 0 0 K ; p c (T=14280 0 K) = 2.3564 x 10' 5 g/cm 3 
(A) Penetration Depth of Keyhole Due to Argon Plasma Jet 
Based on Lancaster (1986) , pressure acting on the free surface of 
the liquid at the bottom of keyhole can be considered such that the 
stagnation pressure due to the plasma jet shall be balanced by the 
surface tension and hydrostatic head, namely 
2a 

P s = + (6-2) 

R 

where P s is the stagnation pressure due to plasma arc jet; p m , metal 
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density of workpiece; <r, the surface tension of liquid metal of 
workpice; and R, the radius of curvature at the root of penetration. 
In the case of VPPA welding system operated at NASA/MSFC, by using AL 
#2219 as a workpiece, values of material properties can be obtained as 
follows (American Society for Metal Handbook Committee, 1979 ; Robert 
and Melvin, 1981) : p m = 2.84 g/cm 3 , and <r (at melting temperature of AL 

#2219 for wall temperature T u = 916. 67°K) = 0.865 N/m. The stagnation 
pressure due to the plasma arc jet can be shown as 



1 (puA) 2 

2 pA 2 



2 P A 2 


Thus, the minimum penetration depth h becomes 

2 o- 1 m 2 2 <s 

R 2 pA 2 R 

h = = (6-3) 

Pm*? Pm 1 ? 


' 1 (9 . 585 x 10* 2 g/s) 2 

.2 (2.3564 x 10" 5 g/cm 3 )( 0 . 12363 cm 2 ) 2 

(2) (0.865 N/m) (10 3 g-m-s' 2 /N)‘ 

(1/2) [ (0.1562) (2 . 54 ) cm] 

/ [(2.84 g/cm 3 ) (9.8 m/s 2 )(10 2 cm/m)]" 1 
= 1.4492 cm = 0.57 inch (6-4) 

In our case, the thickness of workpiece is 0.375 in. which is smaller 
than 0.57 in. In other words, the thrust force induced by the plasma 
arc jet can penetrate the workpiece completely. 

Note that the result of penetration depth shown in this analysis 
is strictly based upon the dynamics of the computation of fluid 
mechanics which can provide the necessary condition, but not the 
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sufficient condition. As the plasma arc temperature is several 
orders of magnitude greater than the melting temperature of metal, it 
can assure that the penetration depth of 0.57 in , which is greater than 
the thickness of workpiece (0.375 in.), can be carried out and the 
plasma arc jet can penetrate the workpiece completely. 

Further discussion of the sufficient condition as to whether the 
power of plasma is large enough to provide the thermal energy for 
melting the workpiece will be accomplished in Section (VI-E) . 

(B) Convective Heat Transfer Loss for Partial 
Penetration of Workpiece 

There are two kinds of convective heat transfer on the workpiece. 
At the beginning before a keyhole is formed, plasma arc jet impinges 
perpendicularly on the surface of the weld pool, and is deflected 
across the weld surface radially away from the stagnation point (see 
Figure 5) . Once a keyhole is formed, a significant amount of plasma 
gas will pass through the keyhole, and discharge to the other side of 
the workpiece. 

(B-l) Weld Pool Convective Heat Transfer 

Governing equation of weld pool convective heat transfer, in 
particular for flow with impingement, can be shown from Metcalfe and 
Quigley ( 1975) , 

Nu = 0.78 Re 0 ’ 5 Pr 0,33 (6-5) 

Again, Nu stands for the Nussett number? Re, the Reynolds number; and 
Pr, the Prandtl number. Radius of the area of impingement is assumed 
to be twice the radius of the plasma arc jet. Average temperature is 
taken between the ethalpys of plasma arc and that of weld pool which 
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are the function of temperature. The enthalpy average is given by 
H(T) = l/2[H arc (T a ) + H veld pool (T wp ) ] (6-6) 

Now, from the available VPPA welding system operated at NASA/MSFC, we 
have 

H a r c ( T a ) = h c = 6858.74 cal/g = 28.716 MJ/kg 
H wp (T wp ) = 0.3672 MJ/kg 

where T wp = 916°K is the melting temperature of workpiece. 

Thus, the enthalpy average is 

H(T) = 1/2(28.716 + 0.3672) = 14.5417 MJ/kg 
From tables shown in Drellishak et al. (1963), we have the average 
plasma arc jet temperature T = 12456°K. Also, we have the values of 
physical parameters at this temperature range as follows: C p = 

5.8461 J/g*°K, p = 33.576 g/m 3 , u = 188.6 x 10 -3 g/m 3 and k = 3.1104 
W/m*K. We can calculate the velocity of plasma arc jet, Reynolds 
number and Prandtl number as follows: 

m (9.585 x 10" 2 g/s)(10" 2 m/cm) 

u = = = 329 m/s 

p g A (2.3564 x 10' 5 g/ cm 3 ) (0 . 12363 cm 2 ) 

puD (33.567 g/m 3 ) (329 m/s) (0 . 1562x2 . 54 cm)(10* 2 m/cm) 

R e " " I 

u (188.6 x 10 3 g/m*s) 

= 232.33 

C p u (5.8461 J/g*k) (188.6 x 10" 3 g/m-s) 

P r = = = 0.35448 

k (3.1104 W/m-K) 

Substituting the values of Reynolds number and Prandtl number in 
the heat transfer equation, we can determine the Nusselt number 
for the heat transfer flow field, 

N u = 0.78R e ° • 5 P r ° • 3 3 = (0.78) (232.33) 0 1 5 (0.35448) 0 • 33 
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8.4142 


From the definition of Nusselt number, the coefficient of heat 
convection can be determined 


h = 


N u k (8.4142) (3.1104 W/m-k)(10' 2 m/cm) 


= 0.65964 W/cm : 


D [ (0.1562) ( 2 . 5 4 ) cm ] 

Thus, the plasma gas power loss for the partial penetration of arc jet 


through the weld pool is 


AH 



= (0.65964 W/cm 2 -K) [V 4 { (0. 1562) (2.54)cm) 2 ] 

• [28.716-0.3672) x 10 3 J/g] • (5.8461 J/g-K) -1 

= 1582 W 


Percentage of plasma arc power loss for partial penetration through 
the weld pool is 
Q 1582 

- = = 23% 

P 6979.13 

(B-2) Convective Heat Transfer Through the Wall of Keyhole 
Recent study indicates that heat transfer rate out of the rear 
half of a keyhole is less than that of the front half. This implies 
that more heat transfer area is required for the front half than that 
of the rear half. In fact, the shape of keyhole with which plasma gas 
actually contact and transfer heat energy is like a finger ring. This 
ring-shaped geometry can be simplified and expanded into the drawing 
shown in Figure 6. 

Surface of keyhole has been equally divided into 5 segments along 
the axis of keyhole between the inlet end and outlet end. Calculation 
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procedures of heat convection for Segment 1 is illustrated. Results 
of the rest of the four segments are shown in the Table 2 . Computation 
methods are all similar to that of Segment 1. 

(1) Geometry for Segment 1 

Figure 6 shows the illustration of geometrical parameters. 
From the available data of VPPA welding system, operated at NASA/MSFC , 
we can compute the geometrical parameters as follows: 


a, = 


( (0.1562) (2.54)] 

= (tt) = 0.6232 cm 


'1 ' 

- l 2 
2 


a 


= tan' 


7t 

- r 

12 ) 


= tan 


- i 


0.375 


\2j 


(0.1562) 


= 56.8 


a 2 = a : + 2- 


= 0.6232 + 2 


tana 0 
= 0.8725 cm 


Area = - (a t + a 2 ) L = 
2 


(0.375) (2.54) 


tan (56.8°) 


(1 


\2 


(0.6232 + 0.875) 


((0.375) (2 . 54) 


= 0.14247 cm 2 

(2) Convective Heat Transfer for Segment 1 

The governing equation of convective heat transfer for 
Segment 1 can be obtained from Metcalfe and Quigley (1975) . It leads 
N u = 0.664 R e ° • 5 P r ° • 33 (6-8) 

where N u stands for average Nusselt number; R e , the Reynolds number; 
and P r> the Prandtl number. 

It is noted that heat transfer equation for plasma flow in a tube 
from Hsu and Rubinsky (1987) was used to compute the convective heat 
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transfer of plasma gas flowing through circular tube-like nozzle, 
shown in Equations (4-3) and (4-8) . There are great differences for 
flow patterns between plasma gas through nozzle and that of plasma gas 
through workpiece. These differences are as follows: (1) Nozzle is 

always a solid circular wall while keyhole of workpiece is constantly 
a wall with molten metal rather than a pure solid wall; and (2) Cross- 
section of keyhole is a curved surface like a finger ring-shaped 
rather than a cylindrical or conical-shaped as that of nozzle. This 
makes that Equation (6-8) is different from Equation (4-3). 
Furthermore, the difference in convective heat transfer equations 
between Equations (6-5) and (6-8) are due to the fact that Equation (6- 
6) is applicable to impingement type flow pattern convective heat 
transfer plasma gas in a molten metal wall of workpiece while that of 
Equation (6-8) is suitable for plasma gas flowing in a molten metal 
wall with completed penetration. 

Plasma gas enthalpy at the inlet end of Segment 1 is equal to the 
plasma gas enthalpy at the entrance of keyhole 
H i,i ~ h c = 28.716 MJ/kg 

Again, the iteration process is needed to carry out by assuming the 
temperature profile of the flow field. Assuming plasma gas enthalpy 
at the outlet end of Segment 1 is 
H t # 0 = 24.1255 MJ/kg 

Thus, the average enthalpy of plasma arc jet between inlet and outlet 
ends of Segment 1 is 

H l , m = 1 / 2 ( H i,o + H i,i) = 1/2(24.1255 + 28.716) 

= 26.421 MJ/kg 
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Furthermore, the average enthalpy of boundary layer between the 
average enthalpy of the main flow of plasma arc jet and the wall 
enthalpy of workpiece is 

H x = 1/2 (H j ^ m + H 1/H ) = 1/2 (26.421 + 0.3672) 

= 13.394 MJ/kg 

The corresponding average temperature for the enthalpy can be 
determined from the table shown in Drellishak et al (1963) which is 
= 12 2 63 0 K 

Physical parameters of flow fields at this temperature are as 
follows: C p = 5.351 J/g- °K, p- 3 4.694 g/m 3 , u = 193.425 x 10" 3 g/m-s, 

k = 2 .937 W/m- °K. 

We can calculate the corresponding Reynolds number and Prandtl number 
as follows: 

mL (0.0985 g/s) [(0.1905 cm)(l0 2 cm/m) 

R e = — = = 76.36 

uA (193.425 x 10“ 3 g/m • s) (0 . 12363 cm 2 ) 

uC p (193.425 x 10“ 3 g/m-s) (5.351 J/g-°K) 

P r = = — = 0.3524 

k (2.937 W/m- °K) 

Substituting the obtained values of R e and P r in the governing 

equation, we have the mean Nusselt number for the flow field 

N u = 0. 664R e ° • 5 P r ° • 33 

= (0.664) (76.36) 0 • 5 (0.3524) 0 ‘ 33 = 4.099 

From the definition of Nusselt number, we can determine the mean value 

of heat covection coefficient 

_ kN u (2.937 W/m- °K) (4 . 099) 

h = = — — ( 10 “ 2 m/cm) 

L (0.1905 cm) 

= 0.632 J/cm 2 -s • °K 
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Plasma arc_power loss at Segment 1 can be determined by 
hA&H 



(0.632J/cm 2 -s-K) (0.14247cm 2 ) [ (26.421 - 0 . 3 672 ) xlO 3 J/g] 


(5.351 J/g- °K) 


= 439 W 

O 

This result is agreeable with the value computed from m(H 1 ^ - H x 0 ) . 
Table 2 shows the computed convective heat transfer loss from Segments 
1 to 5 which are checked with the corresponding enthalpy difference 
between the inlet end and outlet end of each corresponding segment. 
(3) Percentage of Plasma Arc Power Loss Through the Workpiece 
Total plasma arc power loss through the workpiece is the 
summation of power loss for Segments 1 to 5, shown in Table 2. Thus, 
the percentage of plasma arc power loss through the workpiece can be 
determined. It is 


P t + P 2 + P 3 + P 4 + P 5 
P 

439 + 497 + 494 + 350 + 211 1991 

= = = 28.5% 

6979 6979 

(C) Temperature at the Outlet End of Workpiece 
From Table 2, we have the enthalpy for the outlet end of Segment 5 
of the workpiece 

H S/0 = 7.923 MJ/kg 

The corresponding plasma arc temperature at the outlet end of 
workpiece, thus, can be determined from the computed enthalpy based on 
the table shown in Drellishak et al. (1963) 
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T 5 ' 0 = 10870 0 K = T d 

(D) Calculation of the Puddle Width 

Surface of the keyhole has been equally divided into 5 segments 
along the axis of keyhole between the inlet end and outlet end. As the 
temperature gradient along the axial keyhole direction of the 
workpiece is usually negligible at every segment, heat flow during the 
time period of welding can be assumed to be two-dimensional. The 
steady state, two dimensional heat flow during the course of welding 
of thin sheet was first analyzed by Rosenthal (1941) . However, the 
temperature gradient along the axial keyhole direction becomes 
important when we are interested in computing crown width and height 
together with root width and height. In other words, three- 
dimensional heat flow shall be considered when one is particularly 
interested in the relation between crown and root sections of welding. 
In this study, a moving coordinate system is adopted and shown in 
Figure 7 . 

Following assumptions are made for the mathematical 
formulations: (1) No melting and no heat fusion is involved; (2) 
thermal properties of the workpiece are constant; (3) there is no heat 
loss from the surface of workpiece due to the radiation or convective 
heat transfers; and (4) a point heat source is assumed. 

The governing equation for three-dimensional heat conduction is 
given by 


a. 


a z T 3 2 T 


+ 


\ 3x ‘ 


3y' 


3 2 T) 


3Z 2 1 


3T 

- U — = 0 
3x 


(6-9) 


where a (= k/pC p ) is the thermal diffusivity of the workpiece; and U, 
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the traveling speed of welding. This equation is subjected to the 
following three boundary conditions: 


3T 

— = 0 at Z = 0 or Z = W 
3Z 


f 3T' 

Lim -27vr 2 Wk — 
r-» o v 3rj 


= Q 


(6-10) 


T = T 0 as r-»<D 

where r = (x 2 + y 2 + Z 2 ) 1 / 2 is the radial distance; W, the thickness of 
the segment; and T 0 , the temperature of the workpiece before welding. 
General solution of Equation (6-9) subjected to boundary conditions 
(6-10) is a complicated summations of exponential functions. Tsai 
(1982) introduced a semi-empirical analytical solution with the 
consideration of simplification of general solution. The solution 
can be given by 


Q 


'ur' 

2 ~ 

exp 

- 



2 7tkr 


^ 2a) 



In order to estimate the crown and root widths, let us assume that 
the point source at upper half of the convective heat loss through the 
workpiece contritutues to the thermal melting process of crown width 
from solid to molten states, while the point source at lower half of 
that provides the thermal energy similarly for the melting process of 
root width. As we can observe from the facts that both the governing 
equation and boundary conditions are linear, solution of temperature 
distribution based on the contributions made by two point sources, 
point source at the upper half of workpiece and point source at the 
lower half of workpiece, shall be a superposition of the solutions 
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from these two point sources. This 


can be shown as 


T - T 0 = 


exp 


Qi 


ur' 


exp 

- 1 


I 2 " 

2 7rkr x 


k 2a; 



+ 


2irkr 2 



r ur 2 ' 


i 2a ; 


2 1 


( 6 - 12 ) 

where Q l and Q 2 denote heat source intensities of upper heat source and 
lower heat source, respectively. Here r x and r 2 stand for the radial 
distances from the point interests to upper heat source and lower heat 
source, respectively. Figure 8 shows the geometrical coordinate of 
this cross-sectional area of welding puddle. Let us consider to 
solve for crown width d c width an assumption of temperature T to be 
melting temperature of metal T m , we have 


Qi 



dc 1 

2 “ 

~ “ — 0Xp 

2irk(d c /2) 

■ 

1 2a 

2 , 



Q 2 

+ exp 

2 7tk[ (d c /2 ) 2 + (W) 2 ] 1 / 2 

From Table 2, we have the upper half of the convective heat loss 
as 

£(Pj) u = 439 + 497 + 1/2(494) = 1183 W 

and the lower half of the convective heat loss as 
E( p j) u = 1/2(494) + 350 + 211 = 808 W 

Also from Table 2 , we have upper half and lower half heat transfer 
surface area of finger ring-shaped molten metal keyhole as 
£( A j) u = ?(Aj) L = 0.14247 + 0.190 + 1/2(0.2256) = 0.44527 cm 2 


f . 

u ' 

( 

\2cl ) 


(6-13) 


[ (d c /2) 2 + W 2 ] 
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Upper half and lower half circumferential area of keyhole is 


A u = A L = ttDW/2 

= tv[(0.1562 in) (2.54 cm/in) ][ (0. 375 in) (2.54 cm/in)/2] 
= 0.5936 cm 2 

and 

A„ 


Qi = 


J 


?(Pj)u 


and 


0.5936 


0.44527 


(1183) = 1577.1 W 


Q 2 = 


? ( A J ) L 1 


f ( p iU 


0.5936 


(808) = 1077 W 


0.44527 

In our case, workpiece material is AL #2219-T87, we have material 
constants and the grouping of physical parameters as follows: 
p = 2.84 x 10 3 kg/m 3 , C p = 864 J/kg-K, k = 130 W/m>°K 
T m = 916 K, U = 8 in/min = 3.387 x 10" 3 m/s 


130 W/m- °K 


a = 


pC p (2.84 x 10 J kg/m J ) (864 J/kg-K) 
Substituting these values to Equation (6-13) , one has 
(1577.1 W) 


= 5.298 x 10' 5 m 2 /s 


618 K = 


exp 


7T ( 13 0 W/m • k) d c 

(1077 W) 


( 

(3.387 x 10" 3 m/s) d c ' 

2x 

i 

_ 4(5.298 x 10' 5 m 2 /s) _ 

/ 


7t ( 13 0 W/mk) [d c 2 + (4) (9.525 x 10" 3 m) 2 ] 1/2 
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• exp 


3.387 x 10" 3 m/s 


12 


[ (4) (5.298 x 10~ 5 m 2 /s J 


(d c 2 + (4) (9.525 x 10" 3 m) 2 ] 


Solving for d c , we have crown width as 
d c = 8.22 x 10 -3 m = 8.22 mm 
Similarly, we have root width as 

d r = 6.58 x 10 -3 m = 6.58 mm 

(E) Power Consumed in the Full Penetration of Workpiece 
Most of the power absorbed by the workpiece in the welding is used 
in the melting of the workpiece. Assuming the puddle is and circular. 
The volume of puddle melted per unit of time is 

w 


1 



2 

f dr l 

2 

[del 



- 

- 7TW 
3 


Lj 

+ 

12 J 

+ 

' 2 J 


1 2 j 



+ ~ (d c + d r ) U 
2 


7T W 

= — W (d c 2 + d r 2 + d c d r ) + - (d c + d r ) U 
12 2 


12 


(9.525 x 10 “ 3 m)[(8.22 x 10 -3 m) 2 + (6.58 x 10 -3 m) 


+ (8.22 x 10" 3 m) (6.58 x 10" 3 m) ] + - (9.525 x 10" 3 m) 

2 

• [(8.22 x 10" 3 m) + (6.58 x 10“ 3 m)] (3.387 x 10 -3 m/s) 

= 6.50 x 10 -7 m 3 /s 

It is assumed that the puddle temperature is the melting temperature 
of the workpiece. The power absorbed by the workpiece for the full 
penetration of plasma arc is 


• « • 

Q = pvH + pvCp (T ffl - T e ) 

where H is the heat of fusion, and T e is the environmental temperature 
of workpiece. In our case, the workpiece material is AL #2219-T87, 
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the corresponding material constants are as follows: 
p = 2.84 x 10 3 kg/m 3 , H = 387 x 10 3 J/kg, 

C p = 864 j/kg • K, T e = 298 K, T„, = 916 K 
Thus, the power absorbed by the workpiece for full penetration of 
plasma arc is 

Q = P v [H + C p (T m - T e ) ] 

= (2.84 x 10 3 kg/m 3 ) (6.50 x 10 -7 m 3 /s)[(387 x 10 3 J/kg) 
+ (864 J/kg • K) (916 -298k)] 

= 1700 W 

From the flow parameters shown in Table 2 and the calculation of crown 
and root widths, the total absorption of power for plasma arc jet 
between the inlet and outlet ends of workpiece is 1991 W which shall be 
able to provide enough energy for the plasma arc jet to penetrate the 
workpiece completely. In this calculation, it shows that thermal 
power needed to melt workpiece is 1700 W. In other words, the thermal 
energy provided by the plasma arc jet which is absorbed by the 
workpiece is far greater than the energy needed to melt the material 
completely for penetrating through the workpiece. 

Summary of flow zones including flow through nozzle, flow 
between the outlet of nozzle to the workpiece, and the flow through 
keyhole are illustrated in Table 3. 

(VII) Calculation of Crown and Root Heights 
Assuming the shape of the bead reinforcement is in the form of 
parabolic. The coordinate system is chosen to establish in crown 
bead. h c and h r are the crown and root heights, respectively. 
Figure 8 shows the adopted coordinate system for the mathematical 
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formulation. 


The configuration of crown bead can be shown in the following 
equation : 

4h„ 


Z = - 


y + h c 


The cross-sectional area of crown reinforcement is 


A c = 2 


Hc/l ' 


' o 


4h, 


dc 2 


y + h c 


dy = - h c d c 

3 


Similarly, the cross-sectional area of root reinforcement is 
2 
3 


A r = — h r d r 


As the cross-sectional area of welding increment is provided by the 
wire feed, following relation is derived based on the conservation of 
metal volume and the gap of workpiece: 

2 2 * V v 

~ d c h c + “ d r h r = -d v 2 W W (7-1) 

3 3 4 U 

where d w denotes wire feeder diameter; V w , the wire feeding velocity; 
U, the traveling seed of welding; and W g , the gap distance between two 
workpieces. Figure 9 shows the free body diagram and the geometry of 
the coordinate. 

From the figure shown in Figure 9, the cross-sectional areas of 
a, b and c illustrated can be calculated separately as follows: 


' d r/ 2 

4h c . 

1 

r dr ] 


-y 2 + h c 

dy = h c d r h c d c 


' 0 

l d t ! J 

3 

^ d oJ 


A b = d r W 


41 



The slopes at the points B and D are 


tan a 


dz 

dy 


y = d r/2 



tan 3 


dz 
dy 

Y = d r/ 


4h r 

d r 


From free-body diagram, shown in Figure 9, the equilibrium equation of 
forces gives 

1 f d r ) 3 2 

h c d r h c d c — + d r W + - d r h r P g 

3 u e ; 3 

4h_ 4h c fd_ \ ' 

<r - Pd r (7-2) 

dp d c v d c J _ 

(A) Order Estimation of External Pressure Acting on the Crown 
It is found from the experimental observation of the crown and 
root heights of welding that the effective external pressure acting on 
the surface of crown is in the pressure zone between the stagnation 
pressure of plasma arc jet and the hydrostatic pressure of shielding 
gas when the VPPA welding system was traveling through the crown 
during the process of material solidification. This assumption can 
be checked from the mathematical computation of the external pressure 
acting on the crown based on the profiles of crown and root measured 
from the experiment. 

Experimental measurements of the profiles of crown and root, as 
shown in the Figure 8, are as follows: 
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8 mm 


d c =9 mm, d r = 
h c = 1.3 mm , h r = 1.3 mm 

and W = 6.35 mm for the material of AL #2219-T87 used for 
workpiece. 

Rewriting equation (7-2) , we have 
flu 


P = 8 


ld r 


\ 

i 


3 2 


<r — 

h c - - 


h c + W + — h r 

pg 

/ 

3 

vd c j 

3 



(7-3) 


Measurement of specimen is based on vertical welding 
condition. Under this conditon, direction of bead heights and the 
direction of gravitational force are mutually perpendicular. This 
makes that the contribution of gravitational force is not effective on 


the determination of bead heights. With neglecting gravitational 
force, Equation (7-3) becomes 


(hr h c ) 


This equation is true only for puddle in a molten state while the 
experimental measurement is made after the solidification of puddle. 
Shrinkage of crown height will be resulted in a phase transition from 
^oiten metal to solidified metal. Detailed consideration of the 
metal shrinkage during the solidification process will be discussed 
in Section VIl-c of this report. Results from Section VII-C, shown in 
Equation (7-15) , will be used in this study for the modification of 
crown height, h c , after the phase transition from molten metal to 
solidified metal. From Equation (7 — 15) , shrinkage of crown height is 

3 aV 3 

4h c '= - W — = - (6.35) (0.067) = 0.32 mm 

4 V 4 
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Here, W stands the thickness of workpiece; and aV/V, the specific 
volumetric shrinkage. King (1987) shows that the specific 
volumetric shrinkage for Aluminum during the solidification process 
is 6.7%. Thus, the crown height before the solidification process is 
given by 


h c '= h c + Ah c '= 1.3 + 0.32 = 1.62 mm 


Now contribution from root surface tension pressure is 

h r 


1.3 x 10' 3 m 


Psr = 8 


d r 2 


a - ( 8 ) 


(8 x 10' 3 m) 2 

= 140.56 p a 

Contribution from crown surface tension pressure is 


(0.865 N/m) 


h, 


1.62 x 10~ 3 


m 


s c 


-<r = (8) 


= 138.40 P. 


(9 x 10“ 3 m) 


(0.865 N/m) 


Thus, the external pressure acting on the crown is 


P = P - P 

r sr r sc 


= 140.56 - 138.40 = 2.16 P. 


(7-5) 


(7-6) 


In other words, the external pressure acting on the crown is 


P 2.16 

= = 1.5% 

P sr 140.50 

of the root surface tension pressure, and is 


P 2.16 

= = 1.5% 

P sc 138.4 

of the crown surface tension pressure. Based on this estimation, it 
is shown that the external pressure acting on the crown is less than 
1.5-s of the either crown surface tension pressure or root surface 
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tension pressure- Thus, it can be ignored from the crown and root 
heights calculation. 


(B) Crown Height and Root Height Under Normal Gravity 
Rewrite Equation (7-2) with neglecting the pressure term, 



= - d r W 


(7-7) 


Solving for h c from Equations (7-1) and (7-3) , we have 



In our case of VPPA welding system operated at NASA/MSFC, following 
values of physical parameters are adopted: <j = 0.865 N/m, g = 9.81 
m/s 2 , V v = 38 in/min, U = 8 in/min, W = 9.525 x 10~ 3 m, and W g = 3.387 x 
10“ 5 in- Aluminum density at 25°C, P - 2.84 x 10 3 kg/m 3 . Aluminum 
density at molten state before the solidification can be calculated by 
taking specific volumetric shrinkage into consideration, P ' = p (l - 
AV/V) = 2.84 x 10 -3 (1-0. 067) = 2.65 x 10 3 kg/m 3 . 


<r 0.865 N/m 

P g (2.65 x 10 3 kg/m 3 ) (9.81 m/s 2 ) 


3.33 x 10" s m 2 


V w it p \ 2 (38 in/min) 

- d w 2 — = — x 0.0254 m = 9.402 x 10‘ 6 m 2 

4 U [4 (16 j _ (8 in/min) 


d r 6.58 mm /’d r ) 3 

— = = 0.8; d c — = (8.22 x 10' 3 m)(0.8) 3 

d c 8.22 mm ld c / 
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= 4.216 x 10* 3 m 

Thus, we can determine the crown height by substituting the values of 
above mentioned parameters in Equation (7-8) 


h c = - 


1 - 


(12) (3.33 x 10' 5 m 2 ) 
(6.58 x 10" 3 ) 2 


[9.402 x 10" 6 m 2 


(9.525 x 10" 3 m) (3.38 x 10" 5 m) ] - (6.58 x 10" 3 m) • 


(9.525 x 10' 3 m) 


/ (6.58 x 10" 3 m) 


(2 X 

- x 8.22 x 10" 3 m 
3 


+ (8) (3.33 x 10" 5 m 2 ) 


0.8 


8.22 x 10" 3 m 


(7-9) 


[ 8.22 x 10" 3 m (6.58 x 10" 3 m) 

- - (4.216 x 10" 3 m ) j = 0.16 mm 

Similarly, we have the root height 
h r = 1.87 mm 

Width: 8.22 mm 

Height: 0.16 mm 

Width: 6.58 mm 

l Height: 1.87 mm 

under normal gravity conditions. 

(C) Crown and Root Heights Under Vertical Welding Condition 

Under vertical welding condition, direction of bead heights 
and direction of gravitational force are mutually perpendicular. 
This makes the contribution aof gravitational force on the 
determination of bead heights ineffective. 
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Thus, for Crown Bead 


for Root Bead 






Let us consider Equation 



(7-2) with neglecting P and g for vertical welding condition. 
Following equation is obtained: 



Based on Robert and Melvin (1981) for the material of AL #2219 , we 
have relation for surface tension coefficient as the function of 
temperature distribution, 


<r = 865 - 0 . 14 (T - T m ) dyne/cm (7-11) 

where T m denotes the melting temperature of the material. 

Assuming that root surface temperature is same as the 
temperature of crown surface, Equation (7-10) can be simplified to 


h r 


(7-12) 


One can solve Equations (7-1) and (7-12) simultaneously with 
substitution of following data based on VPPA welding system operated 
at NASA/MSFC: d u = 1/16 in, U = 8 in/min, V v = 38 in/min, W g = 3 . 387 x 

10" s m, d c = 8.22 mm, and d r = 6.58 mm. It leads 
h c - 1.09 mm 


and 


h r = 0.70 mm 

These crown and root heights stand for heights of crown and root at the 
stage of molten metal state. There is a phase transition from liquid 
to solid during the metal solidification which accompanies with 
shrinkage of volume. 

Figure 10(A) shows the configuration of the cross-section of 
puddle with the interface between liquid (molten metal) and solid 
metal. As the temperature start to drop down, solidification of 
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metal is observed at the interface between the molten and solid metals 
of the workpiece. Figure 10(B) shows the configuration of puddle 
with a trapezoid-shaped liquid-solid interface which moves gradually 
toward the center due to the solidification of metal. It shows that 
the solidified metal is accumulated toward the root side because root 
temperature is lower than the crown temperature which makes the 
solidification faster in the root side than the crown side. Figure 
10(C) shows the configuration of the cross-section of puddle during 
the process of solidification. This figure shows that the shrinkage 
of crown height is observed during the process of metal 
solidification. Cross-sectional area of puddle with molten metal 
can be approximated as 

1 

A : - Wd e ' 

2 

where W stands the thickness of workpiece; and d c ' , the crown width in 
a state of molten metal. Let us choose specific volumetric 
shrinkage, av/v, to be a parameter of material property during a 
metalic solidification process, variation of cross-sectional area of 
puddle can be shown in the following relation: 


aA = A 


' ZV'' 

1 

' AV' 


= - Wd c ' 


J 

2 

lv ; 


(7-13) 


It is shown in Equation (7-1) for the cross-sectional area of root 
reinforcement as 


A A = A 


- d c 'h' 

13 


2 

3 


= - d c ' 4h' 


(7-14) 


From Equations (7-13) and (7-14), it leads 
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f *v\ 


1 

2 


Wd c ' 


\v 



Ah' 


Thus, 


3 AV 

Ah' » - W — (7-15) 

4 v 

Based on King (1987) , specific volumetric shrinkage for Aluminum 
during the solidification process is 6.7%. Therefore, we have the 
crown height shrinkage 


(3) 


Ah' 


(9.525 mm) (0.067) = 0.48 mm 


u; 


This calculation implies that the final crown height after the metal ic 
solidification becomes 1.09 - 0.48 = 0.61 mm. Thus, we have 

( width = 8.22 mm 1 

crown 

l height = 0.61 mm 


(7-16) 


width = 6.58 mm 


root j 

(height = 0.70 mm / 

This result shows that root height is greater than crown height even 
for the case that we neglect the contribution of gravitational force 
for the vertical welding if we take the shrinkage effect during the 
metal solidification into consideration. 

VIII. Discussion and Conclusion 
Flow profile and power distribution of argon plasma gas as a 
working fluid in the VPPA welding process has been analyzed. The 
purpose of this study is to provide a suitable algorithm for 
computation simulation of the flow field, energy distribution and the 
prediction of the final production of welding in terms of crown and 
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for flow of plasma jet between outlet of nozzle to workpiece is 2 0%; 
and for flow through keyhole of workpiece is 23% with partial 
penetration, and 28.5% with full penetration. Majority of power 
absorbed by the keyhole of workpiece is used for melting the solid 
metal workpiece to become molten metalic puddle. 

For the heat transfer effect on the workpiece, recent study 
imdicates that heat transfer rate out of the rear half of a keyhole is 
less than that of the front half. This implies that more heat 
transfer area is required for the front half than that of the rear 
half. In fact, the shape of keyhole with which plasma gas actually 
contact and transfer heat energy is like a finger ring. 

Penetration depth of keyhole due to argon plasma jet has been 
calculated based on hydrodynamic balance of stagnation pressure of 
plasma jet, surface tension and hydrostatic head, and also the amount 
of thermal power absorbed by the workpiece to be used for melting 
keyhole with the power provided by the arc plasma jet. It shows that 
current VPPA welding system operated by the NASA/MSFC is capable to 
make full penetration for workpiece with a thickness of 0.375 in. 

Crown and root widths, and corwn and root heights have been 
calculated. Crown and root widths are the function of thermal 
diffusivity, traveling speed of wilding, and melting temperature of 
metal. Crown and root heights are the function of wire feeder 
diameter, wire feeding velocity, traveling speed of welding, gap 
distance between two workpiece, crown and root widths, and the method 
of welding , such as vertical or horizontal welding. Configuration 
of crown and root reinforcement is the equilibrium of surface tension 
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pressure, gravitaitonal pressure and external pressure acting on the 
crown. Experimental measurement based on the welding speciment 
received shows that the order of external pressure and gravitational 
pressure acting on the crown is much smaller than that of surface 
tension pressure and gravitational pressure effective to the 
configuration of crown and root reinforcement. Root height is always 
greater than crown height for the case of horizontal welding under 
normal gravity condition. For the case of vertical welding condition 
in which the gravitational force is ineffective on the determination 
of bead heights, it shows that crown height is greater than root height 
when the puddle was in molten state before the phase transition from 
liquid (molten) to solidified metal. Solidification of metal starts 
at the interface between molten and solid metals of workpiece. It 
shows that the solidified metal is accumulated toward the root side 
because root temperature is lower than the crown termperature which 
makes the solidification faster in the root side than in the crown 
side. The volumetric shrinkage, during the phase transition from 
liquid to solid, makes the final figure of root height greater than 
that of crown height after the metal solidif icaton . 

Present research provides an algorithm for promoting automatic 
control of flow parameters and the dimensions of final product of 
welding specification to be used for the VPPA welding system operated 
at NASA/MSFC. This algorithm is essential for the automation of VPPA 
welding system, in particular for the operation of effective welding 
in a hostile environment such as Space Station at low Earth orbit. 
With the regulation of power input, such as pilot and main arc currents 
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and voltages for straight and reverse polarities, volumetric flow 
rates of argon gas for working fluid and helium gas for shielding gas, 
traveling speed of welding, wire feeder diameter, wire feeding 
velocity, traveling speed of welding, etc., one can very well modify 
the end product of welding through the computation of algorithm 
provided in this study. 
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Figure Captions 


Figure 1. 


Figure 2 . 


Figure 3 . 


Figure 4 . 


Figure 5 . 


Figure 6 . 


Figure 7. 


Figure 8 . 


Figure 9 . 


Figure 10. 


VPPA electric power system. Left hand side shows the 
straight polarity mode, and the reverse polarity mode is 
shown on the right hand side. 

Electric power input system and the dimension of the 
geometry for VPPA welding process. 

Plasma arc power distribution in the whole VPPA welding 
system. 

Electric current and voltage drop distribution of the 
VPPA welding system. 

Impingement of plasma arc jet on the surface of weld 
pool . 

Ring-shaped geometry of simplified heat transfer area of 
keyhole and it dimensions. 

Moving coordinate system for three dimensional heat 
flow. X-coordinate is along welding traveling 
direction; Y-coordinate transverse to traveling 
direction; and Z-coordinate along the keyhole 
direction. Welding traveling speed is - U. 
Geometrical coordinate of the cross-sectional area of 
welding puddle. 

Free body diagram and the geometry of coordinate to 
describe the balance of forces among external pressure 
P, surface tension force, <r, and gravitational force g. 
Configuration of the cross-section of puddle with the 
liquid (molten metal) -solid metal. (A) Molten state 
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puddle before phase transition, (B) Molten state puddle 
in the process of phase transition from molten state 
metal to solidified state metal; and (C) solidified 
metal of welding puddle. 
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Table 1: Power Loss Between Plasma, from Outlet End of Nozzle to Plasma 

at the Inlet End of the Keyhole of Workpiece 



Total 













Table 2 Convective Heat Loss Through Workpiece 



Total 1991 1993 


































Table 3 Flow Zone Characteristics 


* 


Flow 

Zone 


Average 

Temperature 


Reynold ! s 
Number 


Prandtl' s 
Number 


Major Loss of 
Heat Transfer 


Percentage of 
Power Loss 


Flow 

Characteristics 


Entrance 

Temperature 


Outlet 

Temperature 


Boundary 

Layer 

Thickness 

Comparison 


Flow Through 
Nozzle 


10,680 °K 


107.6 


0.76 


Convective 
Heat Transfer 


29.6% 


Laminar 

Flow 


12,600 0 


8,760 Q K 


Subsonic 


K 


Thickness 

(Momentum Boundary 
Layer) 

Thickness (Thermal 
Boundary Layer) 


Flow Between 
Outlet of Nozzle 

to Workplprp 


10,720 °K 


139.4 


0.604 


Radiative 
Heat Transfer 


20 % 


Turbulent Subsonic 
Flow * 


8,760 °K 


11,260 °K 


Thickness 
(Momentum Boundary 
Layer) ^ 

Thickness 
(Thermal Boundary | 
Layer) 


Flow Through 
Keyhole 


12,575 °K 


165.7 


0.36 


Convective Heat 
Transfer 


Partial 

Penetration 


23% 


Full 

Penetration 


28.5% 


Laminar Subsonic 
Flow 


14,280 °K 


10,870 °K 


Thickness 

(Momentum Boundary 
Layer) ^ 

Thickness (Thermal 
Boundary Layer) 


Flow becomes Turbulence 


for free 


stream jet with ReNo. exceeding 100 
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